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Abstract. The Landau-Lifshitz-Gilbert equations for the evolution of the magnetization, in 
presence of an external torque, can be cast in the form of the Lorenz equations and, thus, 
can describe chaotic fluctuations. To study quantum effects, we describe the magnetization 
by matrices, that take values in a Lie algebra. The finite dimensionality of the representation 
encodes the quantum fluctuations, while the non-linear nature of the equations can describe 
chaotic fluctuations. We identify a criterion, for the appearance of such non-linear terms. This 
depends on whether an invariant, symmetric tensor of the algebra can vanish or not. This 
proposal is studied in detail for the fundamental representation of u(2) = u(l) x su(2). We find 
a knotted structure for the attractor, a bimodal distribution for the largest Lyapunov exponent 
and that the dynamics takes place within the Cartan subalgebra, that does not contain only 
the identity matrix, thereby can describe the quantum fluctuations. 


1. Introduction 

Recent advances in magnetic materials and techniques allow manipulation of spin moments at 
nanoscale resolution. Thus, chaotic fluctuations become signihcant and their control have been 
the subject of both experimental and theoretical extensive studies [1]. One of the first direct 
observation of period doubling and chaos was spin-wave instabilities in yttrium iron garnet (YIG) 
and has been documented more than thirty years ago [2]. Using the ferromagnetic resonance 
technique, several routes to chaos have been found and explored including periodic-doubling 
cascades, quasi periodic and intermittent dynamics which exhibit complex magnetic behaviors 
[3l Sj. Recently, the phase diagram of a chaotic magnetic nanoparticle has been presented [5], 
which was obtained by monitoring the classical dynamics of its magnetization, modeled by the 
Landau-Lifshitz-Gilbert (LLG) equation. 

What has received much less attention is the contribution of quantum fluctuations, that 
become significant at nanoscale resolution and are crucial for controlling qubit devices [6] . These 
might affect non-linear effects in new ways. 

The challenge, therefore, is to describe the interplay between the two sources of fluctuations in 
a way that can lead to a deeper understanding of their effects and predict new features. To this 
end a model is proposed that displays both and allows us to distinguish them in a particularly 
clean fashion. 


Let take a closer look to the LLG equation. The motion of the magnetization M in 3D space 
that spins around an effective pulsation oj with a damping vector D and an external torque T, 
is given by 


M = u X M - 3 + f. 


( 1 ) 


This equation can be used to describe a particular magnetic system with 3 single-ion anisotropy 
axes n, of intensity r/j, under a static magnetic field B and a diagonal Bloch-Blombergen damping 
D = {Mi/ti, M 2 /T 2 , M^/ts) with the torque T = (0, 0, d) along the z-axis. In such a situation, 
w = XlLi 'ni{ni-M)ni- j3 r]i = (5 = B/ixqMs, Ki is the anisotropy energy along 

the i-axis and Mg is the saturation magnetization. Equation 0 can be written in components 
as 

Ml 

Ml = (?72 — %)M2M3 — /32M3 -I- /I3M2- 

Ivl 

M 2 = ivs —'ni)M3Mi — PsMi + - - (2) 
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M 3 = {r]i - ri 2 )MiM 2 - /3iM2 -|- l32Mi -h d 
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The linear transformation 


Ml = X 

M2 = y 

M 3 = z — r — a 


(3) 


when 7]! =2,r]2 = 1,??3 = = 0,^2 = 0,^3 = cr,n = 1 /(t,T2 = 1,T3 = 1/6, d = -b{r + a), lets 

the Landau-Lifshitz-Gilbert equations Q take the form of the well known prototypical system 
that displays the full repertoire of behaviors from regular to chaotic, along many routes, namely 
the Lorenz system [7] 

X = a{y — x) 

y = x{r - z) -y (4) 

i = xy — bz 


whose typical solution (in the chaotic phase) is displayed in Fig.[^ Here x = dx/dt and similarly 



Figure 1. (Golor online) Plot of the parametric solution of the time evolution of 
Tr(X), Tr(y) andTr(Z) with cr = 10, r = 28, 6 = 8/3. Glassical Lorenz system in gray and 
its matrix version in red. 

for the other variables. While these equations were deduced to describe classical fluid dynamics 
and the original parameters reflect this fact: a is the Prandtl number, r the Rayleigh number 
and 6 is the aspect ratio of the “cell” (real or virtual), their scope is, in fact, much broader, as 


was realized from the work in the 70’s [H]. These are, still, classical equations and the variables 
(x, y, z) are commuting quantities. 

To describe the quantum fluctuations, a generalization of these equations to the case where 
the variables become operators, {X,Y,Z) has been considered in ref. [9], specifically, the model 
discussed in ref. [T^, where the operators are described by square matrices, that are expanded 
in the generators of a given Lie algebra. As is usual in quantum mechanics, a prescription 
for the expressions that involve products of non-commuting variables is mandatory. The Weyl 
ordering m is adopted, so that the product XY of two operators, X and Y is replaced by its 
symmetric expression XY —)• | (XY -|- yA)h The X, Y, Z are now expanded in the generators, 
r“, of a Lie algebra, G, where a = 1, 2,..., dimAd(e), thus X = Y = y°'T'^ and Z = 

A Lie algebra is defined by its structure constants, that enter in the commutation relations, 
[r“,T^] = and the fact that this algebra is compact implies that Tr(T“r^) = k5°‘^ with 

K 7 ^ 0 and Tr stands for the trace on the algebra. 

Here we come upon the following subtlety: The Lie algebra we are considering here, must 
contain a u(l) factor. The reason is that, otherwise, the equations are inconsistent, since the 
trace of the left hand side vanishes, but the trace of the right hand side does not: the trace of 
the non-linear combinations is non-zero. 

Taking into account such requirements, the quantum counterparts to eqs. Q take the 
following, intriguing, form 

(5) 

= -62“ + 

These equations, herald the appearance of the invariant symmetric tensor 



of the Lie algebra. This tensor appears in gauge theories, since the gauge fields belong to the 
adjoint representation of the group and plays an important role in the classification of gauge 
anomalies [12]. Groups, for which this tensor vanishes identically are called “anomaly-safe” 
and in the present context, this means that eqs. ([^ are linear, thus do not give rise to chaotic 
fluctuations. Only quantum fluctuations can appear. Groups, for which this tensor does not 
identically vanish, on the other hand, lead to non-linear equations, thus can describe both 
chaotic and quantum fluctuations. 

For the case of abelian Lie groups, the structure constants vanish. In that case, eqs. Q 
is equivalent to eqs. Q up to a rescaling of all the variables proportional to the single non¬ 
vanishing element of d. This means that they share the same route to chaos, i.e. belong to the 
same universality class. 

When the d tensor no longer vanishes, chaotic fluctuations and quantum effects mix in a non¬ 
trivial way. The simplest case is the u(2) = (u(l) xsu(2))/Z2 algebra, where all the components of 
d“^^ are zero, except d*^“ = 1/2 for c G {0,1, 2, 3} with the corresponding circular permutations. 
This can be made particularly clear by writing the corresponding Lorenz system of equations as 
follows: 

^ It may be checked that another ordering prescription does not change, qualitatively, the results presented here. 
A more complete discussion will be presented elsewhere. 




u(l) 


= a{y^ — x^) 

yO — _yO _|_ _ 2 - 0^0 _ 2 x^z^ 

= —bz^ + + 2x^y^ 


r x“ = a{y°‘ — x“) 
su(2) < = —y“ + rx“ — 

I i;“ = —62;“ + x“y° + x°y“ 


( 7 ) 


The first set of the three equations highlights the fact that u(l) is responsible for the chaotic 
fluctuations and the second set shows that the su(2) components satisfy linear equations, with 
the u(l) variables acting as sources. 

The two sets of equations were integrated using an eighth-order Runge-Kutta scheme with a 
hxed time stepping scheme. The group invariants, Tr(X), Tr(y) and Tr(Z), plotted in Fig. 
define a subspace where a structure similar to the Lorenz “butterfly” appears, “decorated” by 
the 5u(2) terms that make it “knotted”. This last property is, in fact, expected [13], given 
that great circles on the su(2) manifold, the 3-sphere, do have non-zero linking number, the 
well known Hopf invariant [Hj. While knots are stable in the non-chaotic phase, since they 
describe periodic orbits, it will be interesting to investigate in more detail the import of the 
matrix structure of the equations on their properties. 



Figure 2. (Color online) Phase diagram of the largest Lyapunov spectra in the u(2) matrix 
Lorenz system. Details of selected parts of the diagram are inserted to exhibit the Icr error bars. 


A striking feature of the largest Lyapunov exponent, as function of r, is its bimodal 
distribution, along the two group factors, u(l) and su(2), in the non-chaotic phase, r < rent, 
whereas, in the chaotic phase, r > rent, the distribution becomes unimodal, giving the same 
values for each group factor. While a mathematical proof for this result is not available, we 
stress that it provides a consistency check for the reliability of our numerical analysis, since the 
su(2) factor cannot give rise to chaos by itself. The transition to chaos, at r = rent, appears at 
the same value as for the classical Lorenz system, only much more abrupt in the matrix case. 

Another way to characterize the quantum fluctuations is by studying the time evolution of 
the three commutators, Tr([A, T]), Tr([y, Z]) and Tr([A, Z]). If all vanish, this means that all 
three matrices belong to the Cartan subalgebra. Preliminary results show that this does, in fact, 
occur uni and our numerical results seem to confirm it. What is noteworthy in Fig. is that the 
time evolution of the average of the three commutators taken over random initial conditions, 
collapses along the same curve. We shall now eliminate the possibility that the Cartan subalgebra 
















Figure 3. (Color online.) Plot of the average of the three commutators 
Tr([X, y]), Tr([y, Z]), Tr([X, Z]) taken over random initial conditions for r = 15. 


so obtained is proportional only to the identity. This is achieved by monitoring the invariants 
in the su(2) subspace, namely TtX^ — (TVy)^ = and so on. Performing runs 

for different values of r, at fixed b and a, we find that, in the non-chaotic phase, r < Xcrit, 
these quantities seem to converge to a single point, not the origin, whereas in the chaotic phase, 
r > ^’crit) they appear to describe a fuzzy region, of finite volume in phase space. That it does 
not collapse to the origin indicates the persistence of the quantum fluctuations in both phases. 

A natural generalization of our results is towards a quantum Nambu description of spin 
systems, where it has been demonstrated that the Nambu mechanics leads to novel identities 
for extended Lorenz system with dissipation that are not obvious in an Hamiltonian approach 
|15j . Moreover this generalization may also provide insights into the origin of dissipation in such 
magnetic systems that has become a subject of topical research m and lead to new relations 
that are hard to guess from the Hamiltonian viewpoint. 
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